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Abstract. 

We consider the Minimal Supersymmetric Standard Model (MSSM) with varying amounts of non-universality in the soft 
supersymmetry-breaking contributions to the Higgs scalar masses. In addition to the constrained MSSM (CMSSM) in which 
these are universal with the soft supersymmetry-breaking contributions to the squark and slepton masses at the input GUT 
scale, we consider scenarios in which both the Higgs masses are non-universal by the same amount (NUHM1), and scenarios 
in which they are independently non-universal (NUHM2). As the lightest neutralino is a dark matter candidate, we demand 
that the relic density of neutralinos not be in conflict with measurements by WMAP and others, and examine the viable 
regions of parameter space. Prospects for direct detection of neutralino dark matter via elastic scattering in these scenarios are 
discussed. 
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INTRODUCTION 

TeV-scale Supersymmetry (SUSY) is one of the most compelling theories of physics beyond the Standard Model for 
several reasons: it facilitates unification of the gauge couplings, as expected in Grand Unified Theories (GUTs); it 
stabilizes the Higgs vacuum expectation value, offering a solution to the related heirarchy and naturalness problems 
of the Standard Model; and it predicts a light Higgs boson, as expected from electroweak precision measurements. In 
addition, if R-parity is conserved, the lightest supersymmetric particle (LSP) is stable, and, if uncharged, is therefore 
a natural particle candidate for astrophysical cold dark matter. 

Phenomenologically, one of the most-studied versions of the Minimal Supersymmetric Standard Model (MSSM) 
is the constrained MSSM (CMSSM), in which the soft supersymmetry-breaking contributions to the masses of the 
SUSY partners of the quarks and leptons and the Higgs scalars are universal at some GUT input scale, as motivated by 
minimal supergravity. However, while the CMSSM may be the simplest scenario, it is by no means the most plausible 
version of the MSSM. Here, we present the results of recent studies of models with non-universal supersymmetry- 
breaking contributions to the Higgs masses |l|,|2|]. 

The CMSSM is parametrized by four continuous parameters specified at the SUSY GUT scale and a sign choice: 
the universal gaugino mass, Wj n, the universal scalar mass, mo, the universal value for the trilinear couplings, Aq, the 
ratio of the Higgs vacuum expectation values, tan(j3), and the sign of the Higgs mixing parameter, jj,. In this scenario, 
the GUT-scale values of the effective Higgs masses are m\ = m% = mo, and the electroweak vacuum conditions fix |ju | 
and the pseudoscalar Higgs mass, m^. In the NUHM1, the effective Higgs masses are assumed to be universal at the 
GUT scale, though the universality with mo is relaxed. One more input parameter, in addition to those of the CMSSM, 
is required, and may be taken as fi, itia, or the GUT-scale value of the Higgs masses mi — 1112- Similarly, in NUHM2 
models, no relation between m\, m.2, and mo is assumed. This scenario may be parametrized by additional inputs fx 
and tha, or by both GUT-scale effective Higgs masses. 

In each scenario, the renormalization group equations of the MSSM are used to determine the low-scale observables, 
and constraints from colliders and cosmology are applied. We assume that the lightest neutralino is the dark matter 
candidate, calculate its abundance, and determine the prospects for direct detection for models which are phenomeno- 
logically and cosmologically viable, i.e. models not excluded by constraints from colliders, and where the upper limit 
on the relic density of neutralinos is respected. Direct detection experiments derive limits on the dark matter-nucleon 
scattering cross section under the assumption that the local density of cold dark matter comes from a single source, 
that is, one particle species is responsible for the entire dark matter abundance, £Icdm- However, this may not be the 
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FIGURE 1. In the left panel, we show the CMSSM (m^jWin) plane for tanjS = 10, Aq = 0, and > with the corresponding 
neutralino-nucleon elastic scattering cross sections as functions of neutralino mass in the right panel. Constraints are as discussed 
in the text. For further details on this and other figures, please refer to 0]. 

case. If the calculated neutralino abundance is less than the WMAP central value 0], Q.cdm, we assume a secondary 
source of cold dark matter, and rescale the neutralino-nucleon cross sections by a factor Q.%/ Qcdm in order to compare 
our calculated elastic scattering cross sections with limits from direct searches. 

THE CMSSM 

In Fig. Q] we present the neutralino-nucleon elastic scattering cross section as a function of neutralino mass for one 
slice of the CMSSM parameter space. In the left panel, we show the (1W1/2 >"*()) plane for tan j3 = 10, Aq = 0, and ii > 0, 
with the corresponding neutralino-nucleon elastic scattering cross sections for phenomenologically-viable regions in 
the panel on the right. The region excluded because the LSP is a charged ? is shaded brown, and, at large mo, that 
where electroweak symmetry breaking cannot be obtained (resulting in unphysical ii 2 < 0), in pink. The red dot- 
dashed contour corresponds to a Higgs mass of 114 GeV. At lower mw 2 the Higgs boson would be lighter, which 
is excluded by its non-observation at LEP [4]. We also plot a black dashed contour for m x ± = 104 GeV, with the 
region at lower m 1( / 2 a l so excluded by LEP. The green shaded region at very low and iuq is disfavored by the 
measured branching ratio for b — > sy J^t], while the light pink shaded region is favored by the measurement of the 
muon anomalous magnetic moment at the 2-a level [6]. Finally, in the turquoise shaded regions, the relic density 
of neutralinos falls within the WMAP range yfl, For the value tanj3 = 10 presented here, the only cosmologically- 
preferred regions are the coannihilation strip, bordering the T-LSP region, and the focus-point region at large mo, where 
il is small and the LSP is a mixed bino-Higgsino state. Over the bulk of the plane, the relic density of neutralinos 
exceeds the value measured by WMAP by more than 2-a. There are, however, slim strips where £l%h 2 is below the 
WMAP range, which lie between the WMAP-preferred strips and the excluded regions they border. 

Turning to the right panel of Fig. [1] we have separated out the cross section for spin-dependent (SD) scattering 
(red/orange) and the scalar or spin-independent (SI) cross section (blue/green). Here and throughout, we assume the 
current limits on the branching ratios of B — > jiji and b — > sy are respected, as well as the lower limit on the chargino 
mass, and differentiate regions in which the scalar Higgs mass is below the LEP limit (green/orange) and above the 
LEP limit (blue/red). Here, the Higgs mass is below the LEP limit in the coannihilation strip at low rri\ n, corresponding 
to low m x . Both the SI and SD cross sections in the upper right panel contain two distinct regions, that arising from the 
focus-point region at m x < 150 GeV and relatively large a, and a more extended region due to the coannihilation strip. 
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FIGURE 2. Here we present the entire potential range of neutralino-nucleon elastic scattering cross sections as functions of 
neutralino mass for the CMSSM, with 5 < tan/3 < 55, 100 GeV < m l/2 < 2000 GeV, < m < 2000 GeV, and -3my 2 <A Q < 
3mj/2 We consider jx < only for tan/3 < 30. On the left are the cross sections for SD scattering, and on the right are those for SI 
scattering. Also shown in the right panel are upper limits on the SI dark matter scattering cross section as specified in the text. 



In the coannihilation strip, 50 GeV < m x < 400 GeV, where the lower limit on m x is a result of the LEP constraint on 
the chargino mass, and the upper limit on m x corresponds to the end-point of the coannihilation strip for tan/3 = 10. 
At larger values of tan/3, we also find a rapid annihilation funnel in the CMSSM parameter space. 

We also plot in the right panel of Fig.Q~]the limits on the SI cross section from CDMS II $7^ (solid black line) and 
XENON10 [8] (solid pink line), as well as the sensitivities projected for XENON100 |9] (or a similar 100-kg liquid 
noble-gas detector such as LUX, dashed pink line) and SuperCDMS at the Soudan Mine [10] (dashed black line). One 
can see that the SD cross sections are generally larger than SI cross sections, though current upper limits are larger 
than 10~ 2 pb. 

At very low mj/2, CDMS II and XENON10 have definitively excluded some of the region where to/, is below 
the LEP limit. We show this explicitly in the left panel by plotting the reach of current and future direct detection 
experiments in the parameter space. Here and in subsequent parameter space scans, we display contours of scalar 
neutralino-nucleon scattering cross section, scaled by Q. x /Qcdm if necessary, of 5 x 10~ 8 pb (solid green contours) 
and 10~ 9 pb (dashed green contours). A cross section of 5 x 10~ 8 pb is currently excluded by XENON10 for m x = 30 
GeV and by CDMS II for m x = 60 GeV, and will be probed by SuperCDMS for m x up to ~ 1000 GeV. Tonne-scale 
liquid noble-gas detectors such as the proposed XENON IT or a similar detector mass for LUX/ZEP will be sensitive 
to scalar cross sections below 10~ 9 pb for all neutralino masses in the range 10 GeV < m x < a few TeV [9]. These 
detectors will be sensitive to cross sections below 10~ 10 pb over much of the preferred mass range m x ~ (9(100) GeV. 

Of course, we would like to know the range of potential neutralino-nucleon elastic scattering cross sections that are 
possible in the CMSSM. In Fig. [2] we scan over all CMSSM parameters, with SD neutralino-nucleon elastic scattering 
cross sections in the left panel, and SI cross sections in the right panel. One can see that, in both the SI and SD cases, 
there is an upper limit on the cross section as a function of neutralino mass. Focusing on the SI elastic scattering cross 
sections, at low m x < 300 GeV, cross sections generally exceed 10~ 9 pb. The largest cross sections, already excluded 
by CDMS and XENON 10, come primarily from the focus point region at large tan/3. 
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FIGURE 3. The left panel shows an NUHM1 (p.,m ) plane faim 1 / 2 = 500 GeV, tan/3 = 10, andA = 0. The right panel shows 
the corresponding neutralino-nucleon elastic scattering cross sections for phenomenologically-viable models. 



NUHM1 MODELS 



In the NUHM 1 , the parameter space is expanded by one dimension, resulting in additional regions where all constraints 
are satisfied: selectron coahhihilation strips and crossover regions. In the left panel of Fig. [3] we show an example slice 
of NUHM1 parameter space, a (/j.mo) plane. The cosmologically-preferred (turquoise) shaded strips, from large to 
small are: rapid annihilation funnels, where m x w m^/2, rising up at \ jx\ ~ 1000 GeV and arcing outwards (note 
that for ji < the entire funnel is excluded by the branching ratio of b — > sy); selectron coannihilation strips bordering 
the black excluded selectron-LSP regions at low mo; stau coannihilation strips between |/i | ~ 1000 GeV and |/i | w 300 
GeV at low mo; and vertical crossover strips at |ju| « 300 GeV, where the LSP is in a mixed bino-Higgsino state. The 
resulting neutralino-nucleon cross sections are shown in the right panel of Fig. [3] Most cross sections are clustered 
at m x ~ 215 GeV due to the fact that the lightest neutralino is bino-like over most of the plane. The smallest cross 
sections shown come from the funnel or other regions where the relic density is suppressed (and therefore the cross 
section is scaled). The largest cross sections come from the crossover strip, near its intersection with the coannihilation 
strip, and for fi > 0. Between the crossover strips, moving to smaller m x decreases, as does the relic density of 
neutralinos, leading to strips in the right panel of scaled cross sections, with m x decreasing to w 100 GeV. 

In Fig.@] we show a different slice of parameter space, this time the (mA,'«i/2) plane, for mo = 500 GeV, tanj3 = 10, 
pi > 0, and Aq = 0. The triangular allowed region is bounded by the constraint on the branching ratio of b — > sy at 
small niA, the appearance of a slepton LSP at large m^, and the absence of a consistent electroweak vacuum at larger 
m.A and smaller m\ii- The diagonal blue line indicates where m x — ma/2. On each side of this contour there is a 
narrow rapid-annihilation funnel strip where the relic density falls within the WMAP range. The funnel extends only 
to m\/2 ~ 1200 GeV, and therefore m x < 550 GeV for this region. There is another WMAP strip in the focus-point 
region close to the electroweak vacuum boundary, where the LSP is more Higgsino-like. The displayed part of the 
focus-point strip is cut off at tha = 2000 GeV, corresponding to m x < 600 GeV: larger values of m x would be allowed 
if one considered larger values of m^. Apart from the region between this strip and the boundary, and between the 
funnel strips, the relic LSP density would exceed the WMAP range. The limit on the branching ratio of b — > sy is 
important at very low tha and mo, but it is the constraint on the Higgs mass (red dot-dashed curve, roughly horizontal 
at mj/2 = 400 GeV) that places a lower limit on the expected LSP mass for the funnel region of m x ~ 160 GeV. 

The green dashed contour in the left panel of Fig. [4] indicating a future sensitivity to spin-independent dark matter 
scattering at the level of 10~ 9 picobarns, runs mainly through the region where the relic density exceeds the WMAP 
upper limit. SuperCDMS at Soudan would be sensitive to much of the focus-point region shown in the panel on the 
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FIGURE 4. The left panel shows an NUHM1 {m A ,m 1 1 2 ) plane for m = 500 GeV, tan/} = 10, and A = 0. The right panel shows 
the corresponding neutralino-nucleon elastic scattering cross sections for phenomenologically-viable models. 

left, but a more advanced detector would be required if the relic density of neutralinos is below the WMAP range. 
Unfortunately, even if neutralinos make up all the dark matter in the universe, much of the funnel region will remain 
out of reach, even to next-generation direct detection experiments. Points associated with the funnel with cross sections 
larger than 10~ 9 pb fail the Higgs mass constraint. 

The right panel of Fig. |4] displays the scattering cross sections for phenomenologically-viable points in the left 
plane. The fact that scalar cross sections in the funnel region are generally smaller than those in the focus-point region 
is reflected in the lower cutoff on the LSP mass, m x < 550 GeV, seen in the right panel for the points with Osi 
few x 10~ 10 pb. In the focus-point region, the LSP is in a mixed state with m x < 600 GeV, and the cross section may 
be much larger: by two orders of magnitude for the SI cross section, and four orders of magnitude for the SD cross 
section. These large cross sections at large m x have no counterparts in the CMSSM, where mixed states in the focus- 
point region are much lighter. Furthermore, in the CMSSM, the focus point is reached only at large mo, whereas here 
it appears even at fixed mo = 500 GeV. CDMS II and XENON 10 are already beginning to probe he large scalar cross 
sections in the focus-point region, with a few points at very low nt\ n already being excluded, as one can see directly 
in the right panefl We also note that there are points in the right panel with very low cross sections even though they 
have m x < 150 GeV, which also have no counterparts in the CMSSM. These points are associated with the funnel and 
their cross sections are scaled down due to the low relic density in that region. 

A new feature seen in this plot is the near-vanishing of the spin-dependent cross section when m x = 350 GeV. This 
feature is associated with the funnel region where elastic scattering cross sections are substantially lower than those 
of the focus-point region. However, near m x = 350 GeV there is the possibility for a complete cancellation in the SD 
neutralino-nucleon elastic scattering matrix elements when the spin contribution from up quarks cancels that due to 
down and strange quarks. The exact position of the cancellation is sensitive to the values of the spin matrix elements 
adopted, but the existence of the cancellation is robust. Of course, as we present here the cross sections for scattering 
on protons, the cross sections for scattering on neutrons will not exhibit a cancellation in the same place. 



We note these points also have m/, < 1 14 GeV and are within the solid green contour in the left panel of Fig.|4]at very low values of m\j2 an d m A- 
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FIGURE 5. The left panel shows an NUHM2 (m l / 2 ,m ) plane for n = m A = 500 GeV, tan/3 = 10, and A Q = 0. The right panel 
shows the corresponding neutralino-nucleon elastic scattering cross sections for phenomenologically-viable models. 



NUHM2 MODELS 

As already discussed, the NUHM2 has two parameters in addition to those already present in the CMSSM, which may 
be taken as free choices of both the quantities Ma and fl. As the number of parameters is relatively large, a systematic 
survey of the NUHM2 parameter space is quite complex. Here we exhibit one plane of parameter space whose features 
we compare with the CMSSM and NUHM1. For further discussions of direct detection cross sections in the NUHM2, 
see llQj]. 

We display in the left panel of Fig. 0a sample NUHM2 (m{/ 2 ,mo) plane with tan/3 = 10 and fixed niA = 500 GeV 
and ji = 500 GeV, which reveals a few interesting new features. As in the CMSSM, there is a region in the plane 
at large m l ^ 2 an d small mo which is forbidden because the lighter stau is the LSP. Just above this forbidden region, 
as in the CMSSM, there is a stau-coannihilation strip. However, jutting up from this strip at m l ^ 2 ~ 600 GeV and 
~ 950 GeV, there are vertical strips where the relic density of neutralinos falls within the WMAP range. The double 
strips at mi/2 ~ 600 GeV form a rapid-annihilation funnel on either side of the (solid blue) contour where tha — 2m z . 
Again, such funnels appear only at large tan/3 in the CMSSM, but the freedom to choose different values of m^ in 
the NUHM2 permits the appearance of a rapid-annihilation funnel also at the low value tan/3 = 10 shown here. The 
other vertical WMAP strip appears because, as m \j 2 increases relative to pi (which is fixed here), the Higgsino fraction 
of the lightest neutralino increases, which in turn increases the annihilation rate, decreasing the relic density. In this 
case, this results in a crossover region when m^ ~ 950 GeV. At slightly larger m.\ n — 1020 GeV, the lightest and 
next-to-lightest neutralinos are nearly degenerate in mass, and rapid coannihilations (through Z-exchange) result in a 
narrow region with suppressed relic density, and therefore a suppressed scalar cross section. 

These novel regions are clearly visible in the right panel of Fig. The elastic scattering cross sections do not vary 
rapidly as m x increases, until the funnel region at m x ~ 250 GeV is reached. The "V"-shaped suppression in the cross 
section arises from the increasing value of mo as one rises up the funnel. The two sides of the funnel approach each 
other as mo increases, eventually joining together and resulting in a minimum value where the two sides of the funnel 
meet (at a value of mo > 2000 GeV). After this excitement, the cross section continues to rise gradually as one follows 
the coannihilation strip, until the crossover strip is reached at m x ~ 400 GeV. Here the cross section decreases again 
as mo increases, to values even smaller than in the rapid-annihilation funnel, before rising again and finally declining 
towards the end of the coannihilation strip. Note that the entire region to the right of the transition strip is viable, 
albeit with a relic density below the WMAP range. The cross section is therefore reduced due to scaling. Because the 
neutralino is predominantly Higgsino-like here, its mass is given by jj. rather than nii/ 2 , and so we see in the right 
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FIGURE 6. We show the entire potential ranges in the NUHM1 and NUHM2 of the SI neutralino-nucleon cross sections with 
respect to neutralino mass. In both plots, we scan 5 < tanjS < 55, 100 GeV < m,\j2 < 2000 GeV, < rao < 2000 GeV, and 
— 3wi/2 ^A()<3ni l j2 - The GUT-scale values of m\ and m2 are in the range (—2000,2000) GeV, with m\ = mi in the NUHM1. 



panel that the largest values of m x correspond to our choice of fixed pi. 

SUMMARY 

Finally, in Fig. [6] we show the entire potential ranges for the SI neutralino-nucleon elastic scattering cross section as a 
function of neutralino mass in the NUHM1 (left) and NUHM2 (right). Because of the additional freedom in the Higgs 
sector, and since the SI elastic scattering cross sections receive important contributions from Higgs exchange, there is 
more variability in the cross sections in the NUHM1 and NUHM2 than in the CMSSM. This is clear from comparison 
with the corresponding CMSSM panel in Fig. [2] Differences are due primarily to the appearance of additional regions 
of parameter space in NUHM models where the lightest neutralino is in a mixed bino-Higgsino state, such as the 
crossover region. 

In both the NUHM1 and NUHM2, we find significantly larger cross sections at larger m x than would be expected 
in the CMSSM. These cross sections may be probed by the next generation of direct detection experiments for m z as 
large as ~ 700 GeV, while in the CMSSM there is no hope for detection with these instruments for m x > 350 GeV. In 
addition, we find viable points at very low m x with rather low cross sections. If nature has given us a neutralino of this 
character, then we may not directly detect dark matter for some time, but the LHC may find supersymmetric particles 
fairly quickly. In either scenario, we will know that the CMSSM is not an adequate description of nature. We look 
forward to results from both collider experiments and the next generation of dark matter direct detection experiments, 
which may give us some interesting hints about the mechanism of supersymmetry breaking. 
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